AD-AQ97  665 
UNCLASSIFIED 


NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA 

INITIAL  VALIDATION  OF  AN  AMBIENT  AIR  QUALITY  MODEL  FOR  NAVAL  AI — ETC(U) 
DEC  80  T  S  DOUGLAS »  D  W  NET2ER 

NPS67-81-002  NL 


AD  A  09 


LEVEL' 


/y  NPS67-81-,0O2 

NAVAL  POSTGRADUATE  SCHOOL 

■  Monterey,  California 


Approved  for  public  release;  distribution 
unlimited. 


Prepared  for:  Naval  Air  Propulsion  Center 
Trenton,  New  Jersey  09628 


4- 


81  4  13  oti 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  J.  J.  Ekelund  D.  A.  Schrady 

Superintendent  Acting  Provost 


The  work  reported  herein  was  supported  by  the  Naval  Air  Propul¬ 
sion  Center,  Trenton,  New  Jersey  as  part  of  the  Naval  Environ¬ 
mental  Protection  Technology  Program. 


Reproduction  of  all  or  part  of  this  report  is  authorized. 


This  report  was  prepared  by: 

TT'S .  Douglas 
LT.  U.  S.  Navy 


Reviewed  by 


*tst-  ^  TYs.tZz, 

W.  FT  PLaTzER,  Chairman 
Department  of  Aeronautics 


Released  by: 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  page  'VFhAn  Deto  Enfrod) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  REPORT  NUMBER  12.  GOVT  ACCESSION  NO. 

NPS67  -  81  -  002  /tp\/bO  77 

3.  RECIPIENT'S  CATALOG  NUMBER 

4  T|Tl£  end  Subtitle) 

INITIAL  VALIDATION  OF  AN  .AMBIENT  AIR 
QUALITY  MODEL  FOR  NAVAL  AIR  OPERATIONS 

5.  TYRE  OF  REPORT  A  PERIOD  COVERED 

Final  Report 

December  19  80 

8.  PERFORMING  ORG.  REPORT  NUMBER 

7.  Ay  ThO#< r) 

T.  S. Douglas  and  D.  W.  Netzer 

•  .  CONTRACT  OR  GRANT  NUMBERS 

9  performing  organization  name  and  aooress 

Naval  Postgraduate  School 

Monterey,  California  93940 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  *  WORK  UNIT  NUMBERS 

N6237679WR00014 

M  CON  TROLLING  OFFICE  NAME  AND  ADDRESS  , 

12.  REPORT  OATE 

December  1980 

Trenton,  New  Jersey  09628 

11.  NUMBER  OF  PAGES 

126 

r5  MONITORING  AGENCY  name  ft  AOONCSV!/  dltfirtu  from  Controlling  Ottlco)  H.  SECURITY  CLASS,  (ot  thto  rmoort) 

Naval  Postgraduate  School 

Monterey,  CA  93940  Unclassified 


DISTRIBUTION  statement  rot  thl •  R*porf) 

Approved  for  public  release;  distribution  unlimited. 


17  DISTRIBUTION  statement  f Oi  <*•  mbmtrmct  entered  In  Block  30,  It  dllloront  from  Report) 


i*  KEY  #O*0S  <  Continue  on  revoree  9 1 dm  it  neceeeery  md  Identify  by  block  number) 

Air  Pollution 

Aircraft 

Model 


20  ABSTRACT  rContlnue  on  reeeree  ride  It  neceeeery  end  Identity  hr  bl ock  number) 

An  initial  validation  of  an  Air  Quality  Assessment  Model  for 
Naval  air  operations  was  conducted  at  NAS  Miramar,  CA.  A  pre¬ 
viously  developed  model  was  updated  to  appropriately  represent 
1978/79  operations  and  then  evaluated  for  prediction  sensitivity 
to  variations  in  meteorological  and  dispersion  model  parameters. 

A  joint  effort  with  the  Naval  Air  Propulsion  Center,  the  Environ¬ 
mental  Protection  Agency/North rup  Services,  Inc.  and  PMTC,  Pt.  Muj 


DO  1473  eoirioN  o»  '  nov  «j  i*  objouete 

5  N  >  |»2-  M«-  MO  I 


UNCLASSIFIED _ 

JECUPITv  CLASSIFICATION  O*  ’’Mil  RAGE  D«.  *»<•»•<<) 


1 


UNCLASSIFIED 


fuCu  —  T*  CL*HI»C*TlQti  QW  Tmj  ftmtrn  Mmtfl  _ 

was  conducted  to  obtain  detailed  data  over  a  one-week  period.  ■ 

Much  of  the  measured  pollution  concentration  data  (including  all  \ 

■  of  the  suspended  particulate  data)  was  lost  due  to  instrument  and 
I  measurement  technique  problems.  This  severely  limited  the  valida-  : 
i  tion  effort.  Comparison  of  model  predictions  with  the  limited 
|  initial  measured  concentration  data  indicated  that:  (1)  predicted 
J  CO  concentrations  were  in  good  agreement  with  measurement,  (2)  pre-j 
i  dieted  NOX  concentrations  from  aircraft  idle/taxi  operations  were  | 
;  too  low,  and  (3)  predicted  total  hydrocarbons  and  particulate  con-  j 
l  centrations  were  too  high  for  aircraft  idle/taxi  operations  and 
!  too  low  for  environ  sources.  Model  predictions  were  significantly  i 
i  improved  by  increasing  engine  RPM  settings  to  above  idle  for  all  j 
I  modes  normally  specified  as  idle.  Model  validation  efforts  would  j 
!  be  improved  if  one-half  integer  stability  categories  could  be 
!  measured  and  used  in  the  model.  J 
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I.  INTRODUCTION 


In  recent  years  several  mathematical  models  have  been 
developed  to  predict  the  atmospheric  dispersion  of  pollutants 
emitted  from  aircraft-related  activities  at  and  around  airports. 
These  models  have  used  the  steady  state  Gaussian  plume  formu¬ 
lation.  The  Gaussian  formulation  is  used  because  it  is  adapt¬ 
able  to  distances  and  pollutant  travel  times  associated  with 
airports.  An  early  contract  sponsored  by  the  U.3.  Environmen- 

I 

tal  Protection  Agency  (EPA)  resulted  in  a  model  being  developed 
by  the  Northern  Research  and  Engineering  Corporation  (Ref.  1). 
This  model  was  later  modified  by  GSOMET,  Inc.  (Ref.  2),  and 
dealt  specifically  with  civilian  airport  operations.  A  more 
recent  model  has  been  developed  by  Argonne  National  Laboratory 
(ANL)  for  the  USAF  and  was  termed  the  Air  Quality  Assessment 
Model  for  Air  Force  Operations  (AQAM)  (Ref.  3)«  This  computer 
model  was  based  upon  an  earlier  TRW  model,  the  Air  Quality 
Display  Model  (Ref.  1+). 

Each  of  the  models  utilizes  a  method  for  solution  of  dif¬ 
fusion  equations  assuming  Gaussian  dispersion  in  both  the  hori¬ 
zontal  and  vertical  directions.  Gaussian  formulation  in  air 
quality  model  calculations  requires  meteorological  inputs 
including  stability  of  the  atmosphere,  mixing  layer  height, 
and  wind  direction  and  speed.  Detailed  pollution  source  data 
are  also  required.  The  resultant  models  consisted  of  emission 
and  dispersion  programs.  AQAM  included  three  major  parts,  a 
Source  Inventory  model  which  yields  annual  emission  at  an 
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activity  by  source,  a  Short  Term  dispersion  model  which  performs 
hourly-averaged  calculations  using  input  dispersion  parameters 
and  a  Long  Term  dispersion  model.  The  models  predict  average 
steady-state  concentrations  during  the  specified  time  interval 
over  a  specified  grid  surrounding  the  airport. 

Model  verifications  have  to  be  conducted  to  test  the  algo¬ 
rithms  and  plume  dispersion  equations.  Initial  efforts  to 
validate  AQAM  were  begun  by  the  Air  Force  at  Williams  AFB, 
Arizona.  Williams  AFB  was  chosen  because  it  was  a  high  traffic- 
volume,  military  airfield  where  accurate  statistics  would  be 
available.  These  statistics  included  aircraft  type,  mix,  and 
activity  schedules  from  which  emissions  input  data  could  be 
calculated  (Ref.  5)*  The  objectives  of  the  validation  effort 
were  three-fold: 

1.  Collect  a  data  base  of  airport-related  air  quality 
measurements  to  evaluate  the  Air  Force  AQAM  model. 

2.  Determine  the  impact  (if  any)  of  airport-related 
activity  on  local  (5  tern  radius)  air  quality. 

3.  Conduct  a  series  of  special  studies  to  provide  infor¬ 
mation  on  horizontal  and  vertical  dispersion  to 
supplement  any  model  revision  by  ANL  (Ref.  6). 

The  Navy  became  interested  in  the  Argonne  model  capabili¬ 
ties  relative  to  Naval  Air  operations.  Under  sponsorship  of 
the  Naval  Air  Propulsion  Center  (NAPC)  Trenton,  N.J.,  the  Naval 
Postgraduate  School  (NPS),  Monterey,  Ca.,  obtained  copies  of 
both  the  Source  Inventory  and  the  Short  Term  models  of  AQAM 
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for  evaluation  and  adaptation  to  Navy  operations.  Upon  com¬ 
pletion  of  modifications,  a  validation  effort  similar  to  the 
one  at  Williams  AFB  was  planned  at  NAS  Miramar,  California. 

The  Source  Inventory  Program,  as  originally  received  from 
ANL,  computes  annual  emissions  of  three  types  of  sources: 
aircraft,  airbase  (non-aircraft)  and  environment  (off-airbase). 
Each  of  these  types  is  further  reduced  by  geometric  configura¬ 
tion  to  either  a  point,  line  or  area  source.  Data  are  input 
to  the  Source  Inventory  program  relative  to  the  type  and  size 
of  source,  location  of  the  emission  plume  in  three-dimensional 
space  and  the  mass  emission  rate  of  each  pollutant  emitted  by 
the  source.  The  model  input  is  often  comprehensive  and  volum¬ 
inous,  leaving  a  great  margin  for  possible  error.  The  program 
calculates  annual  emissions  and  provides  a  qualitative  ranking 
of  the  contributions  to  the  ambient  air  pollution  of  any  indi¬ 
vidual  source.  It  also  prepares  a  data  bank  containing  source 
characteristics,  annual  emission  rates  and  temporal  distribu¬ 
tion  activity  for  utilization  by  the  Short  Term  program. 

The  Short  Term  program  receives  the  above  compiled  annual 
results  and  calculates  the  dispersion  of  generated  pollutants 
over  a  specified  receptor  grid  during  a  given  hour,  day  and 
month  utilizing  average  meteorological  data  input  for  that 
hour  (Ref.  ?).  For  point  and  area  sources  this  is  accomplished 
by  using  initial  source  dimensions  and  meteorological  stability 
criteria  to  project  a  pseudo-upwind  point  source.  Line  sources 
are  generated  along  the  route  of  travel  of  the  source  vehicles. 


The  ohort  Term  model  utilizes  a  line  dispersion  theory  devel¬ 
oped  by  ANL.  The  line  of  finite  cross-section  is  segmented 
into  shorter  lines,  or  "puffs",  which  are  then  dispersed  from 
pseudo-upwind  line  sources  in  much  the  same  manner  as  point 
and  area  sources  (Ref.  3,3). 

Principal  modifications  to  AQAM  were  required  by  the  Navy- 
due  to  differences  in  flight  operations  between  the  Navy  and 
Air  Force.  Subroutines  were  added  to  AQAM  to  account  for 
Visual  Flight  Rule  (VFR)  approaches  including  aircraft  entry 
break  above  the  runway,  Navy  touch-and-go  cycles,  field  carrier 
landing  practices  (FGLP),  takeoff  delays,  and  hot  refueling 
(refueling  of  aircraft  while  engines  are  operating).  Also 
AQAM  was  expanded  to  handle  helicopter  operations.  It  should 
be  noted  that  modifications  were  only  made  to  subroutines  in¬ 
volving  aircraft  sources.  Airbase  and  environ  source  data 
remain  relatively  consistent  from  base  to  base  whether  Navy, 

Air  Force  or  civilian.  The  Short  Term  portion  of  AQAii  was 
modified  to  calculate  dispersion  of  pollutants  over  h.12  grid 
receptors  rather  than  the  Air  Force's  312  receptors.  This  was 
done  so  that  a  larger  off  airbase  area  could  be  included  in 
the  analysis.  Finally,  Navy  aircraft  engines  and  fuel  types 
are  often  different  than  those  of  the  Air  Force  and,  conse¬ 
quently,  aircraft  performance  data  and  emissions  data  had  to 
be  input  to  reflect  the  changes.  A  plot  routine  was  also 
incorporated  into  AQAH  so  that  predicted  pollutant  distribution 
patterns  could  be  more  readily  observed  (Ref.  9). 


The  aforementioned  model  verification  performed  by  the 
Air  Force  at  Williams  AFB  involved  13  months  of  continuous 
air  monitoring  during  the  period  June  1976  through  June  1977. 
Air  quality  data  were  collected  at  five  ground  stations  and 
meteorological  data  were  taken  routinely  at  the  base  weather 
station.  Aircraft  operations  data  and  airbase  and  environ 
source  information  were  then  input  to  A«iAM  and  predicted 
values  of  pollutant  concentrations  were  compared  with  observed 
or  measured  data  from  the  monitoring  stations.  Preliminary 
results  have  indicated  that  a  reasonable  correlation  exists 
between  predicted  and  observed  hourly  pollutant  levels  (Ref. 
10). 

The  Air  Force  effort  included  a  wide  range  of  meteoro¬ 
logical  conditions  collected  over  a  long  period  of  time.  It 
wa3  decided  to  concentrate  the  Navy  validation  effort  on  a 
specific  meteorological  "window"  which  would  be  reasonably 
stable  for  several  days  and  which  would  occur  when  a  large 
amount  of  aircraft  activity  occurred.  The  latter  was  neces¬ 
sary  in  order  to  minimize  the  problem  inherent  with  high  back¬ 
ground  pollution  levels.  Specifically,  it  was  desired  to 
perform  the  validation  effort  at  NAS  Miramar,  CA  and  to  obtain 
more  detailed  data  relative  to  (1  )  aircraft  taxi  and  refueling 
operations,  (2)  hourly  aircraft  flight  activity,  and  (3) 
meteorology. 

Once  the  Navy  modifications  were  completed  and  input  data 
were  obtained  for  NAS  Miramar,  it  was  necessary  to  determine 
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the  sensitivity  of  the  model  predictions  to  the  input  meteo¬ 
rological  and  operational  conditions  and  to  certain  dispersion 
model  parameters.  An  initial  model  sensitivity  study  was 
performed  using  the  Navy  version  of  AQAM  and  1975  activity  at 
NAS  Miramar  as  a  representative  data  base  (Ref.  9). 

The  purposes  of  the  present  study  were  (1 )  to  update  the 
data  in  the  Source  Inventory  program  of  AQAM  in  order  to  rep¬ 
resent  1978/1979  operations  at  NAS  Miramar  and  (2)  to  compare 
the  predicted  and  measured  levels  of  pollutant  concentrations 
for  the  purpose  of  validating  the  Short  Term  program  of  AQAM. 

A  necessary  component  of  the  validation  effort  was  the  con¬ 
ducting  of  an  updated  model  prediction  sensitivity  study. 

II.  OVERALL  MODEL  VALIDATION  EFFORT 


The  Navy  version  of  the  AQAM  model  validation  effort  was 
initiated  by  the  Naval  Air  Propulsion  Center  (NAPC).  NAPC 
provided  the  funding  and  necessary  program  coordination  as 
well  as  technical  assistance  in  selection  of  the  monitoring 
site  locations  and  the  required  data  acquisition.  NAS  Miramar 
was  chosen  because  it  had  the  largest  number  of  flight  opera¬ 
tions  of  any  NAS  and  because  it  had  been  used  in  previous  work 
performed  by  the  Naval  Postgraduate  School  in  developing  the 
Navy  version  of  AQAM. 

The  overall  objectives  of  the  NAPC  program  were  to: 

a.  validate  the  AQAM  model, 

b.  document  the  effects  of  aircraft  operations  on 
air  quality,  and 
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c.  assess  the  possibility  of  using  AQAM  (as  an  alter¬ 
native  to  an  expensive  monitoring  program)  to 
determine  the  effects  of  aircraft  operations  on 
air  quality  at  other  NASs  (Ref.  11). 

The  program  was  divided  into  two  related  parts.  The  first  part 
was  to  consist  of  a  one  year  continuous  monitoring  study.  In 
late  1979  -  early  1980  air  quality  was  being  measured  24  hours 
a  day  using  an  automated  data  acquisition  system.  This  effort 
was  directed  primarily  at  objective  (b)  noted  above.  The  second 
part  was  to  consist  of  two  special  studies,  each  one-week  in 
duration.  The  latter  studies  were  to  be  intensive  in  nature 
with  detailed  operational,  meteorological  and  pollution  concen¬ 
tration  data  being  collected.  These  studies  were  to  be  directed 
primarily  at  objectives  (a)  and  (c)  above.  The  first  special 
study  took  place  in  August  1979  and  data  received  from  that  week 
were  used  in  the  model  validation  discussed  herein.  The  second 
special  study  was  scheduled  for  the  spring  of  1980  but  was  sub¬ 
sequently  cancelled.  The  two  periods  were  chosen  to  occur  during 
distinctly  different  meteorological  conditions,  especially  lid 
height  and  stability  category.  Organizations  involved  in  the 
special  study  and  individual  responsibilities  of  each  included: 

a.  Northrup  Services  Incorporated  (MSI)  contracted 
by  EPA:  Air  quality  monitoring  and  data  reduc¬ 
tion  to  provide  hourly  averaged  pollutant  levels. 

b.  Pacific  Missile  Test  Center  (PMTC) :  Meteorological 
measurements  and  data  reduction  to  provide  hourly 


averaged  weather  conditions  throughout  the 
receptor  grid. 

c.  NAPC/NPS:  Aircraft  activity  monitoring. 

d.  NPS:  Reduction  of  aircraft  activity  data  for  in¬ 
put  into  AQAM,  model  predictions  using  items  b. 
and  c.  above,  and  comparison  of  predictions  with 
measured  values. 

III.  NAS  MIRAMAR  INTENSIVE  DATA  ACQUISITION 

Planning  the  special  study  for  validation  of  AQAM  began 
with  identifying  both  the  emittants  to  be  monitored  to  best 
characterize  dispersion  and,  as  previously  mentioned,  locating 
appropriate  monitoring  stations. 

The  major  pollutants  in  aincraft  engine  exhausts  include 
particulates/smoke  (PT),  carbon  monoxide  (GO),  unburned  hydro¬ 
carbons  (HG)  and  oxides  of  nitrogen  (NOX).  The  relative  amounts 
emitted  depend  primarily  upon  the  engine  thrust  setting.  In 
addition,  sulfur  oxide  emissions  (SOX)  are  often  significant 
from  industrial  and  domestic  furnaces.  Therefore,  CO,  HC,  NOX, 
PT  and  SOX  were  selected  as  the  pollutants  to  characterize 
emissions  of  both  aircraft  and  airbase  related  activity.  Figure 
1  identifies  the  grid  system  used  to  locate  the  receptors  in 
AQAM.  The  grid  spacing  was  1  km.  and  the  x-y  coordinates 
varied  from  (0,0)  to  (24,1 5)  representing  400  separate  recep¬ 
tor  locations. 


S 


IFICATION  AT  NAS  MIRAMAH 


Locating  measuring  stations  where  continuous-air-monitor- 
ing  instruments  would  be  placed  was  of  prime  importance  in  the 
validation  effort.  The  behavior  of  the  model  predictions  at 
a  particular  receptor  will  depend  to  a  great  extent  on  its 
location  relative  to  numerous  sources  throughout  the  receptor 
grid,  especially  those  located  upwind.  To  validate  the  model, 
it  was  important  to  compare  air  quality  samples  at  locations 
where  the  airbase  and  aircraft  contributions  were  large  rela¬ 
tive  to  background  levels  of  pollution.  Ultimate  placement  of 
the  stations  assumed  a  dominant  wind  from  the  WNW  (292°)  as 
advised  by  PMT C . 

Up  to  12  special  receptor  locations  can  be  input  to  the 
Short  Term  program.  Special  receptor  locations  were  assigned 
to  each  of  the  four  pollution  monitoring  stations  as  indicated 
in  Table  I.  They  are  also  identified  in  Figure  1 . 


TABLE  I 

MONITORING  STATION  LOCATIONS 


TRAILER  NUMBER 

GRID  COORDINATES 

SPECIAL  RECEPTOR 
NUM3ER  IN  A 4AM 

1 

10.01  , 

9. 24 

401 

2 

10.52, 

8.46 

402 

3 

11.24, 

8.35 

406 

4 

12.82, 

7.31 

410 

The  intended  use  of  trailer  1  was  to  determine  background 
levels  of  pollution  upwind  of  aircraft/airbase  sources. 
Trailer  2  was  located  just  downwind  of  the  hot  refueling  site. 
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Trailer  3  was  situated  just  upwind  of  the  hot  refueling  pits. 

It  was  also  downwind  of  the  hot  refueling  area.  Trailer  1; 
was  located  well  downwind  at  the  outer  boundary  of  NAS  Miramar. 

During  the  planning  stage,  N3I  made  equipment  preparations 
for  each  trailer  site  for  the  air  quality  monitoring  experi¬ 
ment.  FMTC  analyzed  the  meteorological  history  for  the  San 
Diego  area  to  determine  the  best  time  period  for  the  special 
study.  Optimum  weather  conditions  for  validation  were  con¬ 
sidered  to  consist  of  a  moderate  wind  coming  from  the  290 
degree  direction,  a  Turner  stability  category  of  2-3»  and  a 
moderate  lid  height  (mixing  layer  depth)  of  lj.00-500  meters. 

It  wa3  desirable  to  have  relatively  constant  weather  condi¬ 
tions  for  the  week  of  intensive  data  acquisition.  This  would 
allow  the  dispersion  model  to  be  validated  with  multiple  tests 
in  which  aircraft  operations  varied  but  weather  remained  approx 
imately  fixed.  The  week  of  1-7  August  1979  was  chosen  as  the 
most  feasible  for  meeting  these  objectives  for  the  first  in¬ 
tensive  study. 

Operating  procedures  for  the  week  proceeded  on  a  previously 
planned  routine.  Specific  tasks  performed  by  NSI  (pollution 
monitoring)  and  PMTG  (meteorological  monitoring)  will  be  pre¬ 
sented  by  those  activities  under  separate  cover.  NP3  and  NAPC 
personnel  monitored  the  detailed  aircraft  activity  in  accor¬ 
dance  with  the  time  schedule  listed  in  Table  II. 


TABLE  II 


AIRCRAFT  ACTIVITY  MONITORING  TIMES  (LOCAL) 


1 

AUG 

1  300 

• 

1600 

2 

AUG 

1000 

• 

1230 

IliOO 

- 

1700 

3 

AUG 

0800 

- 

1230 

6 

AUG 

0900 

• 

1230 

1330 

- 

1630 

7 

AUG 

0830 

• 

1100 

1330 

- 

1630 

Observation  of  aircraft  activity  was  performed/recorded 
from  three  locations  —  the  control  tower,  the  hot  refueling 
site  (octagon)  and  the  refueling  pits. 

The  functions  performed  in  the  control  tower  involved 
(1 )  timing  the  sequences  of  every  aircraft  on  departure  from 
initial  startup  to  takeoff,  (2)  timing  the  sequences  of  every 
aircraft  on  recovery  from  entry  into  the  airport  traffic  area 
(defined  here  as  having  a  three-mile  radius)  to  landing  and 
taxi  to  the  refueling  area.  Also,  the  parking  areas  and  taxi' 
way3  used  by  each  aircraft  and  the  type  of  landing  performed 
(VFR,  IFR)  were  monitored.  Data  sheets  used  to  record  the 
aircraft  activities  observed  from  the  control  tower  are  pre¬ 
sented  in  figures  2  and  3. 

Data  collected  at  the  hot  refueling  sites  (octagon)  and 


refueling  pits  included  time- in-mode,  amount  of  fuel  taken, 
and  aircraft  type  (see  data  sheets  in  figures  U  and  5). 


TAKEOFF  DATA  SHEET 

DUTY  RUNWAY _  WIND _  TIMS 

register/time 

Side  number _ Aircraft  type _ 

Parking  area 

Commence  sequence  0/0 

Start  complete  / 

Taxi  complete _ / 

(holding  at  runway) 

(engine  check  complete) 

Takeoff  complete  / 

EVOLUTION  (check  one) 

Takeoff  and  depart  area 
FCL? 

Touch  and  go 


number 

number 


FIGURE  2 
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LANDING  DATA  SHEET 
( full  atop  landings  onl y ) 

DUTY  RUNWAY  _  WIND _  TIME 


register/time 

Side  number  Aircraft  type 

Commence  sequence _ 0/0 

(enter  break  or  3  mi.  on  IPR  approach) 

Landing  complete _ _ 

(clear  of  runway! 

Taxi  complete 

(pits/hot  refuel  holding  area) 

Fuel  commence _ _ 

(enter  pits/hot  refuel  area) 

Fuel  complete _ _ 

(depart  pits/hot  refuel  area) 

Shut  down _ 

Parking  area _ 

(hot  refuel  aircraft  only) 


FIGURE  3 
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HOT  REFITS L  S3CTJBNCB  DATA  SHBBT 
(OCTAGON) 


A/C  Type  (circle  one) 


Side  number _ 

Arrival  time 
at 

holding  area _ 

Arrival  time 
into 
oc  t  agon 

Departure  time 
from 
octagon 

Pound 3  fuel 

received _ 

Fuel  spilled  yea/no 


F-k 

P-5 

F-S 

?-1  k 

A-k 

3-2 

Transient 


FIGURE  k 


IS 


TIME 


PIT  REFUEL  SEQUENCE  DATA  SHEET 


A/C  Type  (circle  one) 


Side  number_ 

Arrival  time 
at 

holding  area 

Arrival  time 
into 

refuel  pit _ 


F-4. 

P-5 

P-6 

?-i4 

a-4 

E-2 

Transient 


Shutdown 

(circle  one) 

Pounds  fuel 
received 

Fuel  spilled  yes/no 


Hot  refuel 


Departure  time 
from 

refuel  pit _ 

Pounds  fuel 

received _ 

Fuel  spilled  yes/no 


FIGURE  5 
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The  aircraf t/airbase  operational  data  that  were  collected 
were  used  as  input  to  the  Source  Inventory  program.  Air 
quality  measurements  (by  NSI)  and  meteorological  data  (by 
PMTC)  were  also  being  collected  during  the  entire  period  of 
observation. 

IV.  AQAM  MODIFICATIONS  AND  SENSITIVITY  STUDY 

A.  MODEL  MODIFICATIONS 

In  order  to  perform  a  mouel  validation,  the  data  input  to 
the  Source  Inventory  program  must  reflect,  as  closely  as  pos¬ 
sible,  conditions  and  emittant  sources  as  they  exi3t  at  the 
time  of  validation.  Therefore,  one  of  the  purposes  of  this 
study  was  to  update  the  data  in  AQAM  to  represent  1978/1979 
operations  at  NAS  Miramar. 

Changes  made  to  the  input  routines  of  the  A QAM  program 
included  data  input  on  the  3-2  aircraft  —  an  addition  at  NAS 
Miramar  3ince  1975*  Parking  area  coordinates,  taxiway  usage 
and  aircraft  landing  and  take-off  operational  cycle  time-in¬ 
mode  (LTO)  were  all  modified  to  accept  3-2  aircraft  activity. 
All  data  were  input  in  accordance  with  guidelines  stipulated 
in  Refs.  7-9  and  12.  Averaged  meteorological  data  were 
changed  to  reflect  1978  figures.  The  annual  amount  of  air¬ 
craft  activity  for  1973,  including  arrivals,  departures,  touch' 
and-go  cycles,  and  FCLF ' s  was  entered  according  to  aircraft 
type.  The  specific  parking  areas  and  taxiways  used  by  each 
aircraft  were  modified.  Other  emissions  information 


(specifically;  fuel  spillage,  training  fires,  environ  land 
use  area  factors,  and  off  base  vehicle  miles  per  year)  was 
either  added  or  updated.  Airbase,  non-aircraft  activity 
modifications  included  changes  in  test  cell  and  run-up  stand 
usage. 

3.  SENSITIVITY  STUDY  PARAMETERS  AND  PREDICTIONS 

With  the  update  completed  and  reflecting  conditions  as 
they  existed  at  the  time  of  the  first  intensive  study,  an 
investigation  was  performed  to  determine  the  sensitivity  of 
the  model  predictions  to  meteorological  and  operational  con¬ 
ditions  anticipated  for  1-7  August  1979  (special  study).  Sen¬ 
sitivity  results  indicate  under  what  conditions  and  at  what 
receptor  locations  the  model  can  best  be  validated.  In  addi¬ 
tion,  these  results  are  needed  before  conclusions  can  be  drawn 
from  the  comparison  of  measured  and  predicted  pollution  levels. 
In  a  model  validation  effort,  predicted  concentrations  are 
compared  to  measured  values  at  specific  receptor  locations. 

When  making  these  comparisons  it  is  necessary  to  know  how  sen¬ 
sitive  the  model  predictions  are  to  the  uncertainties  in  the 
specified  meteorological  and  operational  input  data.  For 
example,  stability  category  is  normally  specified  as  an  inte¬ 
ger  value  between  one  and  six;  if  the  hourly  averaged  value 
can  only  be  specified  as  two  or  three,  what  effect  would  this 
variation  have  on  the  model  predictions?  In  addition,  it  is 
necessary  to  know  whether  the  monitoring  stations  are  located 
in  regions  where  there  are  large  horizontal  gradients  in  pol¬ 
lution  concentrations. 
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Twelve  special  receptors  were  used  to  examine  the  sensi¬ 
tivity  of  predicted  pollution  levels  in  the  vicinity  of  the 
four  monitoring  stali'ns  to  various  meteorological  conditions 
and  model  parameters.  A  previous  model  sensitivity  study  had 
been  conducted  by  Netzer  (Ref.  9)  using  1975  operational  data 
and  different  nominal  meteorological  conditions.  Table  III 
describes  the  special  receptor  locations  used  in  AOAM  for  both 
the  sensitivity  study  and  the  validation  effort.  Locations 
relative  to  runways,  taxiway3  and  refueling  areas  are  depicted 
in  Figure  1 . 


TABLE  III 

SPECIAL  RECEPTOR  LOCATIONS 

AQAM  RECEPTOR  NUMBER  DESCRIPTION/LOCATION 


u.01 

trailer  1 

k02 

trailer  2 

U03 

100  m  downwind  i 
trailer  2 

of 

kOk 

100  m  crosswind 
of  trailer  2 

( south) 

k05 

100  m  crosswind  (south¬ 
east)  of  trailer  2 

koe 

trailer  3 

k07 

100  m  downwind  of 
trailer  3 

k06 

100  m  crosswind 
of  trailer  3 

( south) 

U09 

approach  end  of 

runway  1 

hi  0 

trailer  If. 

411 

500  m  upwind  of 

k 

trailer 

ki2 

100  m  crosswind 
of  trailer  k 

(north) 
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In  order  to  perform  the  sensitivity  study  it  was  neces¬ 
sary  to  establish  a  reference  or  nominal  case  meteorologically 
and  operationally.  The  anticipated  weather  conditions  for  the 


intensive  study  period,  listed  in  Section  III,  were  used  as 
the  reference  weather.  Meteorological  parameters  were  varied 
independently,  with  aircraft  activity  kept  constant.  Table 
IV  indicates  the  meteorology  data  input  for  each  of  nine  com¬ 
puter  runs. 


TABLE  IV 

METEOROLOGY  FOR  SENSITIVITY  STUDY 


Run  Turner  Wind  Wind  Temperature  Lid 

Number  Stability  Speed  Direction  Height 

(m/s)  (deg)  1  "  (m) 
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Run  number  1  was  the  reference  case.  The  ambient  air  temp¬ 
erature  was  not  varied  because  previous  results  (Ref.  9)  had 
shown  it  to  have  little  effect  on  predicted  pollution  levels. 

The  aircraft  activity  data  input  to  the  Source  Inventory 
program  were  representative  of  one  hour  of  daytime  flight 
operations.  In  addition,  airbase  and  environ  sources  were 
kept  constant  with  updated  19?8  data.  In  the  normal  mode  of 
utilization  of  A QAM,  annual  totals  are  input  and  frequency 
factors  are  used  to  determine  the  total  operations  in  any  one 
month,  week,  day,  and  hour.  For  the  present  effort,  the 
"desired"  one  hour  input  data  had  to  be  scaled  up  to  annual 
operations  in  order  that  the  Short  Term  and  Source  Inventory 
programs  would  function  properly.  The  "scale-up"  factor  used 
was : 

12  hr/day  x  31  day/mo  (Aug)  x  12  mo/yr  =  klx6L  nr/yr  (1  ) 

(12  hr/day  represents  no  night  operations) 

Table  V  presents  the  aircraft  activity  values  which  were  held 
constant  for  the  entire  sensitivity  study. 


TABLE  V 


AIRCRAFT  ACTIVITY  FOR  SENSITIVITY  STUDY 


1  HOUR  OPERATIONS 


AIRCRAFT 

ARRIVALS 

DEPARTURES 

TOUCH 

& 

GO'S 

VFR 

ARRIVALS 

FCLP'S 

F-4 

3 

3 

2 

2 

6 

F-6 

1 

1 

1 

1 

0 

E-2 

1 

1 

1 

1 

0 

F-II4. 

3 

3 

2 

2 

6 

A-1+ 

2 

2 

1 

1 

0 

F-5 

1 

1 

0 

1 

0 

TRANSIENT 

1 

1 

0 

0 

0 

H-3 

0 

0 

0 

0 

0 

1  YEAR  OPERATIONS 


AIRCRAFT 

ARRIVALS 

DEPARTURES 

TOUCH 

& 

GO'S 

VFR 

ARRIVALS 

FCLP'S 

F-4 

1  3392 

1  3392 

3928 

3923 

26784 

F-8 

1+1+61+ 

U14-614. 

W+&1+ 

I+I16I1 

0 

E-2 

1+1+61+ 

1+1+61+ 

I+I4.6I4. 

1+1+61+ 

0 

F-1  k 

13392 

13392 

8928 

8928 

26784 

A-li 

3928 

8928 

I4I4-6I4. 

1+1+61+ 

0 

P-5 

UJU.6I4. 

I4I4-6I4- 

0 

1+1+61+ 

0 

TRANSIENT 

1+1+61+ 

1+1+61+ 

0 

0 

0 

H-3 

0 

0 

0 

0 

0 

Aa  explained  in  Section  I,  the  results  from  the  Source  In¬ 
ventory  program  are  used  along  with  the  meteorological  data  as 
input  to  the  Short  Term  program.  Output  from  the  Short  Term 
program  wa3  arranged  in  3even  tables.  Four  tables  consisted 
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of  pollutant  levels  in  micrograms  per  cubic  meter  from  environ 
airbase,  aircraft  and  total  sources  at  all  specified  grid  re¬ 


ceptors.  2ach  table  listed,  for  all  receptors,  the  receptor 
number  and  x-y  coordinate  location,  and  the  concentrations  for 
all  five  pollutants.  The  remaining  three  tables  expressed  the 
3ame  results  in  terms  of  fractions  of  the  total  emissions  from 
environ,  airbase  and  aircraft  sources. 

The  receptors  of  interest  in  the  sensitivity  study  were 
the  twelve  special  receptors  (i4.OI-l4.i2)  and  that  one  where  the 
maximum  concentrations  existed. 

To  compare  the  expected  affects  of  the  meteorological 
variables  cn  the  predicted  ground  level  (z=0)  concentrations, 
the  Gaussian  dispersion  formula  for  point  sources  can  be  used 
(Ref.  13). 


where : 


*3 

X  =  concentration,  g/mJ 

Q,  =  uniform  emission  rate,  g/sec 

<r  ,  a  a  standard  deviations  of  plume  concentrations 
in  the  horizontal  and  vertical  directions 
respectively,  m 

U  =  mean  wind  speed,  m/sec 

H  =  initial  plume  height,  m 

y  =  0  along  plume  centerline 


V/hen  vertical  diffusion  is  limited  by  a  stable  layer  at 
height  h^^  the  diffusion  equation  is  modified  as  follows: 


For  infinite  line  sources  Turner  (Ref.  13)  utilized: 


X(x,y*z=°;H) 


(U) 


where : 

q  =  source  strength  per  unit  distance,  g/sec-m 

♦  =  angle  between  line  source  and  wind  direction, 
!i5°<  $  <  90° 


Major  variations  of  the  Short  Terra  program  predictions  under 
different  meteorological  conditions  should  follow  equations 
(2),  (3),  or  (ij.) ,  depending  upon  the  receptor  location  rela¬ 
tive  to  the  dominant  emission  sources  (Ref.  9). 


G.  EFFECT  OF  METEOROLOGICAL  PARAMETERS  ON  MAXIMUM  RECEPTOR 

CONCENTRATIONS 

Table  VI  presents  the  predicted  maximum  concentrations 
of  four  of  the  five  pollutants  and  the  location  of  each  for 
the  reference  case.  Also  shown  are  the  maximum  predicted  CO 
and  PT  from  aircraft  sources  for  each  of  the  other  conditions 
investigated.  The  meteorological  variable  is  listed  in  each 
case. 

The  reference  case  indicated  that  the  maximum  contribu¬ 
tions  from  the  environ  sources  occurred  south  of  the  airbase 
(at  receptors  (9,2)  and  (11,2)),  However,  high  levels  of  environ 
pollution  (background)  also  were  predicted  to  occur  throughout 


the  airbase.  On  the  airbase,  the  contribution  from  airbase 
sources  was  generally  negligible,  whereas  aircraft  sources  of 
PT  were  dominant.  Maximum  concentrations  from  aircraft  sources 
occurred  for  CO  and  PT  at  receptor  (11,3),  near  the  intersection 
of  the  runways.  This  was  generally  the  case  for  all  the  con¬ 
ditions  investigated, 

1 •  Stability  Category 

Increasing  the  stability  category  (more  stable  condi¬ 
tions)  decreases  o'  and  ir  ,  and  therefore  should  increase  the 

y  z 

predicted  ground  level  concentration  along  the  wind  vector 
downwind  of  the  source  (see  equation  (2)).  At  the  peak  concen¬ 
tration  receptors  (Table  VI),  which  are  necessarily  near  the 
plume  centerline,  the  increase  in  stability  category  increased 
the  concentration  and  shifted  the  maximum  concentration  receptor 
downwind. 

2.  Lid  Height 

As  a  plume  develops  downwind  of  a  source  it  will  spread 
in  a  vertical,  a3  well  as  horizontal,  direction.  The  ground 
and  lid  height  (elevated  inversion  layer)  act  as  reflectors  of 
the  plume.  Increasing  the  lid  height  would  decrease  the  con¬ 
centration  only  at  receptors  which  are  far  enough  downwind  from 
the  source  for  reflections  to  occur  (see  equation  (3)).  ?or 
the  maximum  receptor  location  (11,6),  lid  height  had  negligible 
effect  on  the  predicted  aircraft  concentrations  of  CO  and  ?T  (Table  VI) 
since  it  was  located  near  the  major  aircraft  sources.  However  the  aircraft 
contributions  increased  with  lid  height  since  concentrations  from  environs 
decreased. 
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Increasing  the  wind  speed  should  decrease  predicted 
concentrations  along  the  plane  centerline  for  a  single  source 
(equations  (2),  (3)  and  (It)).  This  behavior  was  apparent  for 
the  maximum  concentration  receptors  (Taole  71,  run  nos.  6,  1 
and  7 ) . 

L.  Wind  Direction 

Changing  wind  direction  changes  the  orientation  of  tne 
plume  dispersion.  As  a  result,  the  receptor  where  concentra¬ 
tions  were  a  maximum  from  aircraft  sources  was  predicted  to 
shift  to  receptor  (10,8)  when  the  wind  direction  became  270° 
(Table  VI,  run  no.  8). 

J.  EFFECT  CP  METEOROLOGICAL  PARAMETERS  ON  CONCENTRATIONS 

AT  SPECIAL  RECEPTORS 

Short  Term  output  for  each  of  the  nine  sensitivity  runs  is 
presented  in  Appendix  A  for  the  special  receptors.  The  refer¬ 
ence  case  (run  no.  1 )  output  includes  receptor  concentrations 
for  environ,  airbase,  aircraft,  and  total  sources  in  yUgn/m'  as 
well  as  fractional  values  for  aircraft  sources.  Receptor  con¬ 
centrations  for  aircraft  sources  (run  nos.  2-~)  are  included 
in  yugm/'n^  and  fraction  of  total.  In  order  to  visualise  varia¬ 
tions  in  pollutant  concentration,  the  overall  grid  system  was 
mapped  with  contour  levels  for  the  sensitivity  study  in  Appen¬ 
dix  3.  Contours  for  the  reference  case  are  included  for  CC 
and  ?T  concentrations  from  airbase,  aircraft,  and  total  sources 
Contours  for  run  nos.  2-q  are  included  for  CC  and  PT  concentra¬ 
tions  from  aircraft  sources. 


Tables  VIIa-d  summarize  the  special  receptor  concentra¬ 
tions  of  CO  and  PT  for  each  of  the  nine  sensitivity  runs.  In 
general,  the  comments  relating  to  the  maximum  receptor  concen¬ 
trations  pertain  to  the  special  receptor  concentrations.  Prom 
a  modeling  standpoint  special  receptor  lj.01  (trailer  1  )  proved 
to  be  well  located  for  the  purpose  of  measuring  background 
pollutants.  A3  can  be  seen  in  Tables  VIIa-d,  very  little  CO 
and  PT  due  to  aircraft  exist  at  receptor  iiOl  .  When  finite 
values  did  occur  (run  nos.  2,  3,  7,  8  and  9)  they  resulted 
from  the  aforementioned  method  of  projecting  area  sources  (in 
this  case  —  the  hot  refueling  area)  upwind  to  pseudo -point 
sources. 

1 .  Stability  Category 

An  increase  in  stability  category  increases  the  down¬ 
wind  concentration  along  the  plume  centerline  from  a  single 
source  since  the  plume  spreads  more  slowly.  Table  Vila  indi¬ 
cates  that  the  area  around  trailer  1  (receptors  ii.02-Li.05) 
receives  emittants  from  multiple  sources  since  the  concentra¬ 
tions  of  CO  and  PT  first  decreased  and  then  increased  with  in¬ 
creasing  stability  category.  These  receptors  are  also  located 
very  near  large  sources. 

CO  and  PT  concentrations  around  trailer  3  (receptors  L06- 
UOd )  were  significantly  higher  than  those  around  trailer  2  due 
to  the  effect  of  an  increased  number  of  plumes  overlapping 
downwind.  Some  multiple/near  source  effects  were  also  evident 
at  this  location.  The  receptor  concentrations  around  trailer 
k  (receptors  ij.IO-l4.i2)  changed  only  slightly  with  variations 
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in  meteorological  conditions  due  to  the  large  downwind  distance 
from  the  primary  sources.  Concentrations  at  receptor  hC9  were 
high  as  expected  due  to  its  close  proximity  to  runway  and  taxi¬ 
way  line  sources. 

2.  Lid  Height 

At  trailers  2  and  3  lid  height  had  no  effect  (Table 
Vllb).  This  was  expected  since  these  locations  are  very  near 
the  sources  of  pollution.  At  trailer  l|.,  which  is  far  down¬ 
wind,  increasing  lid  height  decreased  concentrations. 

3.  Wind  Speed 

As  indicated  in  Table  Vile,  an  increase  in  wind  speed 
decreased  the  concentration  downwind  at  trailers  3  and  ii. 
igain,  however,  at  trailer  2  the  behavior  was  more  random. 

It.  Wind  Direction 

Changing  the  wind  direction  from  the  reference  290° 
to  310°  (run  no.  9)  resulted  in  the  expected  reduction  in  air¬ 
craft  CO  and  ?T  at  trailers  2  and  3  (Table  Vlld).  In  this  case 
the  plumes  from  the  major  upwind  aircraft  sources  miss  receptors 
ij.02  and  Il06.  However,  when  the  wind  direction  was  changed  to 
270°  (run  no.  6),  the  concentrations  increased  significantly. 
This  indicates  that  trailer  2  was  apparently  outside  the  plume 
from  the  hot  refueling  area  when  the  wind  was  from  290°. 

Further  evidence  of  this  was  that  receptors  I4.OL.  and  I4.05  (cross- 
wind  to  I4.C2 )  had  significantly  higher  concentrations  than  re¬ 
ceptors  kC2  and  L03. 

The  trailer  k  receptor  exhibited  an  increase  in  concentra- 


icn  with  an  increase  in  wind  direction 


This  was  expected 
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since  moat  aircraft  source  plumes  (including  the  maximum  re¬ 
ceptor  location  at  coordinate  (11,3)  are  located  upwind  of 
trailer  lj.,  from  the  290°-310°  direction. 

5.  Special  Receptor  Locations 

As  discussed  above,  for  model  validation  efforts  it  i3 
necessary  to  know  whether  the  monitoring  stations  are  located 
in  regions  where  there  are  large  horizontal  gradients  in  pol¬ 
lution  concentration  or  where  the  concentrations  are  very 
sensitive  to  the  specified  hourly-averaged  meteorological  con¬ 
ditions.  Table  VIII  presents  a  summary  of  the  effects  of 
distance  from  the  monitoring  stations  on  the  predicted  pollu¬ 
tion  concentrations.  Concentrations  are  presented  for  each  of 
the  nine  cases  for  conditions  100m  downwind  and  1 OCm  crosswind. 
As  a  receptor  is  moved  toward  a  specific  plume  centerline,  the 
concentration  would  increase.  When  a  receptor  is  located  down¬ 
wind  from  several  sources,  horizontal  movement  of  the  receptor 
may  increase  or  decrease  the  receptor  pollution  level,  depending 
on  the  multiple  plume  effects. 

Increases  in  concentration  varied  by  factors  of  two  to 
sixteen  at  trailers  2  and  3  for  the  reference  case  as  a  result 
of  moving  the  receptor  100m  downwind  or  closer  to  plume  center- 
line.  No  appreciable  horizontal  gradients  in  concentration 
existed  around  trailer  I4..  In  almost  every  case  (variation  of 
meteorological  parameters),  concentrations  increased  as  expected, 
since  the  receptors  were  moved  closer  to  the  centerlines  of  the 
major  aircraft-related  plumes  for  the  290°  wind.  In  run  no.  8, 
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TABLE  VIII 

DIFFERENCE  FACTORS  IN  SPECIAL  RECEPTOR  CONCENTRATIONS 


Trailer  2 


Trai  ler 


Trailer  4 


lOOir 


100m 


100m 


100m 


HUN 

down- 

cross- 

down 

cross- 

NO. 

wind 

wind 

wind 

wind 

403/402 

404 

402 

40  7 

40C 

103 

406 

i 

CO 

me  10 

inc 

16 

inc 

3 

inc 

2.3 

Xo 

change 

Reference 

pT* 

inc  2  . 

inc 

8 

inc 

4. 3 

inc 

O 

2 

CO 

inc  1.5 

inc 

1 . 5 

inc 

2.5 

inc 

2  3 

Stability 

Category 

PT 

r.o  change 

inc 

2.5 

inc 

3.3 

me 

3T3 

.No 

change 

3 

CO 

me  1.5 

inc 

1.5 

inc 

1.3 

i  no 

1 . 5 

PT 

dec  1 . 1 

inc 

1.  7 

inc 

2.5 

inc 

«">  T 

w  .  xj 

4 

CO 

inc  10 

inc 

16 

inc 

3 

me 

2.  3 

Lid  Height 

5 

PT 

CO 

inc  2 
inc  10 

me 

me 

8 

16 

inc 

inc 

4.3 

3 

inc 

inc 

3 

2. 3 

Xo 

chance 

PT 

inc  2 

inc 

3 

inc 

4 . 3 

i  no 

3 

6 

CO 

inc  1.3 

inc 

1 . 7 

inc 

3.  3 

me 

2.  S 

Wind  Speed 

PT 

CO 

dec  1.5 
inc  1.3 

inc 

inc 

1.8 

1.8 

inc 

inc 

5 

2.5 

i  nc 
inc 

4 

2 

Xo 

on an ge 

PT 

inc  1.5 

inc 

4 . 5 

1  no 

3.  8 

i  ac 

2 . 8 

8 

CO 

dec  1.3 

inc 

5 

inc 

3.5 

i  nc 

1 . 3 

Wind 

Direction 

PT 

dec  2.3 

inc 

2.5 

inc 

4 . 5 

inc 

2 

No 

change 

9 

CO 

inc  2.5 

inc 

1 

inc 

IS 

inc 

13 

PT 

inc  4 

inc 

30 

me 

/  .  O 

i  n  c 

O 

I 
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where  the  wind  direction  was  changed  from  290°  to  270°,  the 
concentration  at  the  100m  downwind  location  decreased  at 
trailer  2. 

These  results  again  indicate  that  comparison  between 
measurements  and  predictions  will  be  most  difficult  at  trailer 

2.  Not  only  do  multiple  plume  effects  and  the  close  proximity 
to  ground  aircraft  sources  cause  unusual  variations  in  con¬ 
centration  with  changing  meteorology  but  also  the  horizontal 
gradients  are  quite  large. 

3.  EFFECT  OF  SPECIFIED  AREA  SOURCE  SIZE  ON  RECEPTOR 

CONCENTRATIONS 

When  large  sources  are  input  into  AQAM  they  are  normally 
modeled  as  area  sources.  The  dimensions  of  the  area  sources 
have  to  be  specified  and  some  judgement  is  required  to  pick 
the  most  representative  dimensions  of  these  "uniform  concen¬ 
tration  sources."  To  determine  what  effect  the  specified 
size  of  aircraft  area  sources  had  on  concentrations  at  various 
receptors,  the  lengths  of  the  sides  of  three  prime  sources 
were  both  increased  and  decreased  by  forty  percent.  The 
specific  sources  included  the  hot  refueling  area,  the  hot  re¬ 
fueling  delay  area  and  the  pit  refueling  delay  area.  The  length 
of  the  sides  of  each  area  source  in  the  reference  case  was  500 
meters.  This  length  was  changed  to  300  meters  and  then  to  700 
meters. 

Increasing  the  size  of  an  area  source  effectively  moves 
the  pseudo-upwind  point  source  further  upwind.  Keeping  the 
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emittants  and  meteorology  constant,  the  plume  would  spread 
at  the  same  rate.  At  a  specific  receptor,  the  concentration 
can  increase  or  decrease,  depending  on  its  location  relative 
to  the  area  sources.  For  this  study,  the  variations  in  con¬ 
centrations  at  trailers  2,  3,  and  k  never  exceeded  six  percent. 

?.  VARIATION  OF  JET  PENETRATION  LENGTH  AND  HORIZONTAL  AND 

VERTICAL  DISPERSION  PARAMETERS 

In  AQAM,  turbojet  exhausts  during  taxi  and  takeoff  are 
treated  as  finite  line  sources.  Initial  line  source  dimen¬ 
sions  and  locations  have  to  be  specified  and  these  are  some¬ 
what  arbitrary.  Currently  in  AZAM  the  jet  is  assumed  to 
"penetrate  1  (i0  meters"  (i.e.,  approximately  U4.O  jet  diameters) 
before  coming  to  rest  relative  to  the  ambient  air.  Default 
values  for  the  line  source  cross-sectional  3ize  are  3m  high 
by  20m  wide.  No  plume  rise  is  considered  to  occur.  These 
line  sources  are  then  treated  as  pseudo-upwind  lines  which 
disperse  in  a  Gaussian  manner  with  the  same  empirical  disper¬ 
sion  oarameters  (o'  ,cT  )  as  used  for  elevated  Doint  sources. 

y  z 

In  a  recent  study  at  the  Naval  Postgraduate  School  (Ref. 

111.)  jet  characteristics  were  measured  in  a  simulated,  neutrally 
stable  atmosphere.  It  was  found  that  jet  penetration  length 
was  considerably  less  than  1UQ  jet  diameters;  being  more  nearly 
35  jet  diameters.  Initial  plume  dimensions  were  found  to  vary 
significantly  with  jet  orientation  to  the  ambient  wind  direc¬ 
tion  and  some  plume  rise  was  observed.  Jet  dispersion  rates 
were  also  found  to  be  more  rapid  than  currently  used  in  AQAM. 


In  order  to  determine  whether  the  above  findings  would 
have  any  significant  effects  on  the  predicted  concentrations 
from  aircraft  sources,  AQAM  was  modified  in  sequential  steps 
as  follows: 

(1)  decrease  the  jet  penetration  length  from  140 
to  35  meters. 

(2)  step  (1)  and  specification  of  initial  aircraft 
line  source  (taxiway  and  runway)  dimensions  as 
a  function  of  orientation  to  the  wind  (per  fig. 
40,  Ref.  14). 

(3)  steps  (1)  and  (2)  and  decrease  the  stability 
category  by  one  to  increase  the  jet  plume 
spreading  rate. 

Decreasing  the  penetration  length  was  found  to  have  little 
effect.  This  was  somewhat  expected  since  the  aircraft  line 
sources  at  NAS  Miramar  have  lengths  up  to  3.7  km.  The  reduc¬ 
tion  in  jet  penetration  length  was  but  three  percent  of  the 
longest  line  source.  In  step  (2)  the  angle  of  incidence 
formed  by  the  wind  with  each  line  source  was  determined,  and 
using  the  a  and  a,  versus  angle  of  incidence  relationship 
determined  by  Brendmoen  and  Netzer,  new  horizontal  and  vertical 
dispersion  parameters  were  input  to  the  Short  Term  program. 

In  general,  the  changes  involved  increases  in  initial  line 
source  dimensions.  At  the  maximum  concentration  receptor  and 
at  trailers  3  and  4,  a  nominal  reduction  in  concentrations  of 
up  to  a  maximum  of  16  percent  was  predicted. 


In  step  (3)  the  above  changes  were  kept  in  AQAM  and  the 
stability  category  was  decreased  from  3  to  2  (more  unstable 
conditions).  Output  indicated  a  decrease  in  concentration  of 
up  to  a  factor  of  two  at  the  maximum  concentration  receptor 
and  at  trailers  3  and  I4..  It  should  be  noted  that  in  its  pre¬ 
sent  form  AQAM  only  allows  variation  of  stability  category  for 
all  dispersions  as  opposed  to  variation  of  aircraft  sources 
alone.  This  decrease  was  expected  as  previously  determined  in 
the  meteorological  sensitivity  study. 

G.  CONCLUSIONS 

Stability  category  and  wind  speed  were  the  two  meteoro¬ 
logical  parameters  that  most  affected  maximum  receptor  concen¬ 
trations.  Model  predictions  will  therefore  be  most  sensitive 
to  uncertainties  in  the  hourly-averaged  values  of  these  para¬ 
meters  which  are  input  into  A QAM.  Wind  direction  had  a  large 
effect  on  the  concentrations  at  trailer  2.  Trailer  2  is 
apparently  located  in  an  area  where  large  horizontal  gradients 
of  pollutant  concentrations  exist,  i.e.,  near  the  edges  of  the 
plumes  from  large  aircraft  sources. 

Trailer  1  appears  to  be  a  good  location  for  measurement  of 
background  pollution  levels. 

Variations  in  aircraft  area  source  sizes  did  not  appreciably 
affect  concentration  levels  at  specific  receptors. 

Variations  of  the  specified  jet  penetration  length  and 
initial  horizontal  and  vertical  dispersion  parameters  of  air¬ 
craft  exhaust  plumes  during  taxi,  takeoff  and  landing  modes 


changed  concentrations  by  a  maximum  of  only  16  percent.  The 
data  of  Brendmoen  and  Netzer  (Ref.  li;)  indicated  that  turbo¬ 
jet  exhausts  spread  more  rapidly  than  elevated  point  sources. 
This  result,  when  incorporated  into  A QAM,  significantly 
affected  predicted  concentration  levels  (by  a  factor  of  2 )  at 
the  monitoring  trailer  locations. 

V.  COMPARISON  CP  AQAM  PREDICTIONS 

WITH 

DATA  FROM  THE  INTENSIVE  STUDY 
A.  VALIDATION  REQUIREMENTS 

As  previously  stated,  model  validation  consists  of  compar¬ 
ing  predicted  hourly-averaged  pollutant  concentrations  to 
hourly-averaged  measured  values  at  specific  receptor  locations. 

A  determination  of  model  accuracy  must  be  made  within  the  con¬ 
text  of  the  accuracy  of  the  input  operational  data  and  of  the 
hourly-averaged  meteorology  and  measured  concentrations.  It 
is  important  to  note  that  although  the  meteorology  and  pollu¬ 
tant  concentrations  may  be  constantly  varying,  only  hourly- 
averaged  values  are  used.  Comparisons  between  measured  and 
predicted  concentration  values  in  areas  where  large  horizontal 
gradients  exist  (trailer  2)  are  likely  to  exhibit  widely-varying 
results.  Because  of  these  factors,  a  need  exists  for  a  vast 
amount  of  accurate  data  with  which  to  conduct  model  validation. 

Prior  to  the  comparison  of  measured  and  predicted  values. 


background  levels/local  air  quality  must  be  determined  in  order 


to  be  able  to  separate  the  contributions  of  aircraft,  airbase 
and  environ  sources  throughout  the  receptor  grid.  The  Source 
Inventory  program  allows  for  input  of  environ  sources.  If 
these  data  are  not  available,  approximate  inputs  can  be  in¬ 
cluded  through  the  use  of  land-use  factors.  The  factors  (Ref. 
12)  distinguish  between  city  center,  urban,  rural,  park  areas, 
etc.  Input  for  off-base  line  sources  (roadways)  requires 
appropriate  vehicle  mileage  and  speed  values.  The  selection 
of  appropriate  land-use  factors  used  in  this  study  was  some¬ 
what  judgemental.  The  roadway  line  source  values  used  were 
based  on  actual  average  daily  traffic  volumes  for  1978  as  pro¬ 
vided  by  the  Comprehensive  Planning  Organization  of  the  San 
Diego  Region.  One  method  for  determining  actual  concentra¬ 
tions  from  aircraft/airbase  sources  is  to  subtract  values  from 
an  upwind  measurement  (i.e.,  trailer  1  data)  from  values 
obtained  at  each  of  the  other  special  receptors. 

Comparison  of  weekend  measured  data  at  each  special  recep¬ 
tor  with  weekday  data  should  also  provide  a  good  indication  of 
background/environ  pollutant  levels  due  to  the  reduction  in 
aircraft  activity  at  NAS  Miramar  on  weekends.  The  measured 
data  indicated  a  wind  speed  varying  from  calm  to  five  knots 
on  Saturday  and  Sunday  approximately  90%  of  the  time.  The 
wind  direction  also  varied  up  to  1 30°  throughout  the  two-day 
period.  This  slight-to-stagnant  air  motion  apparently  caused 
an  accumulation  of  pollutants  at  NAS  Miramar  from  environ 
(local  San  Diego)  sources.  Unfortunately,  this  behavior 
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invalidated  anv  comparison  between  weekday  and  weekend  concen¬ 
trations  for  the  purposes  of  validating  weekday  background  levels 
on  the  airbase.  Therefore,  a  need  exists  for  additional  weekend 
data  when  the  meteorological  conditions  are  more  representative 
of  those  experienced  during  the  period  of  intensive  measurement. 

3.  DATA  REDUCTION'  AND  MODEL  INPUTS 

Measured  data  for  CO,  NOX  and  THC  were  provided  by  NSI  in 
parts  per  million  (ppm).  Comparison  of  these  values  to  AQAM 
predictions  requires  conversion  to  micrograms  per  cubic  meter 
(^gm/m~>'l  .  An  accurate  conversion  exists  for  CO  under  standard 
conditions;  1111.11  x  ppm  CO  =yygm/m'3  CO.  The  most  often  used 
conversion  for  NOX  is  based  upon  NO-,:  2000  x  ppm  NOX  =  xgm/m^ 

NOX.  Measured  data  were  obtained  for  THC  and  CH, .  CH,  often  is 
the  major  portion  of  THC  concentration  in  urban  atmospheres  of 
North  American  latitudes.  Typical  concentrations  are  1.2S-1.5 
ppm  (Ref.  6).  The  CH^  conversion  is  666.6?  x  ppm  CH^  =^gm/m'> 

CH_j .  However,  aircraft  generated  hydrocarbons  may  be  significant¬ 
ly  heavier  than  CH^.  The  only  PT  data  available  were  measured 
by  a  nephelometer  in  terms  of  the  scattering  coefficient,  b  (bscat) 
Air  samples  were  also  taken  to  determine  total  particulates  (TP) , 
but  the  data  were  invalidated  as  a  result  of  a  filter  preparation 
error  by  contractors  at  U.  C.  Davis.  This  loss  of  the  TP  data 
severely  affected  the  model  validation  effort.  Particulate  con¬ 
centrations  on  the  airbase  had  the  best  chance  to  be  dominated 
by  aircraft  sources  and  therefore  provided  one  of  the  best  means 
for  comparing  predictions  with  measurements.  For  the  bscat  data, 
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•in  average  c.mverf  i  an  tact  or  was  emp  1  Jvei  R  e  : I  '  .  '  •  .  . 

N'eph  ;':’scat  1  ^  ^  gm-  mJ  PT .  "he  latter  :or.v?rs :  n  fa;:  jr  ta'.  "e 
ir.  considerable  error. 

The  AQAM  model  was  run  over  ten  one  -  hour  time  periods  as 
listed  in  Table  IX.  The  type  of  aircraft  activity  varied  con¬ 
siderable  throughout  the  ten  VOAM  runs.  ’/.'hen  .1:  f  ferent  from 
normal  operations,  remarks  of  the  act :  vi  *  v  3  ra  included  m  Tab  1  e 
IX.  The  chosen  periods  of  time  were  primaril'.  in  the  at  t  err.  on 
when  the  wind  speed  and  lid  height  are  greatest. 

Figures  6a-e  present  the  meteoro logica  1  conditions  at  NAS 
Miramar  (obtained  from  NAS,  Pt .  Mugu  investigators’  for  the  iav 
of  intensive  measurements  and  detailed  observation  of  aircraft 
activity.  The  values  are  hour  1  -  average i  and  plotted  over  f  e 
1000-2000  time  period  for  each  day .  .Ml  weather  conditions  were 
averaged  over  the  applicable  time  periods  shown  cross  -  hatched 
in  Figures  6a-e. 

Runs  5  and  9  were  performed  to  determine  whether  or  net 
transit  time  of  emittants  affected  predicted  concent  rat  1 'n  level 
relative  to  runs  4  and  S.  A  fifteen  minute  emittant  travel  time 
was  chosen  due  to  the  wind  speed  and  average  distance  from  four; 
to  monitoring  station.  It  should  be  noted  that  the  final  runs 
(6-7  August)  had  significantly  higher  wind  speed,  temperature  m 
lid  height.  This  variation  in  meteorology  was  not  anticipate! 
and  was  somewhat  undesirable  from  a  model  validation  viewpoint. 

Due  to  variations  in  meteorology  within  the  calculated 
dispersion  times,  it  is  generally  agreed  that  the  values  of 
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I'POKOLOGICAL  CONDITIONS  AT  NAS  MI  HAMAH  VS.  TIME  OP  DA 


METEOROLOGICAL  CONDITIONS  AT  NAG  MIRAMAR  VS.  TIME  OF  DAY  (2  AUG) 
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CT  and  d,  cannot  be  more  accurate  than  a  factor  of  2.  In 
addition  to  this  uncertainty,  model  predictions  are  sensitive 
to  the  average  meteorology  used  as  input  as  discussed  above. 

For  example,  consider  the  data  presented  in  Fig.  6a  for  the 
period  1300-1400  hrs.  During  this  period  the  wind  direction 
changed  from  270°  to  300°  and  the  stability  category  changed 
from  1  to  2.  The  lid  height  and  wind  speed  were  steady.  Values 
employed  for  wind  direction  and  stability  category  for  this  period 
(Table  IX,  run  no.  1)  were  290°  and  2,  respectively.  The  sensi¬ 
tivity  study  of  section  IV  has  shown  that  a  decrease  in  wind  direc 
tion  of  20°  and  a  decrease  in  stability  category  from  2  to  1  can 
increase  the  predicted  concentrations  at  trailer  3  by  factors  of 
1.3  and  1.3  respectively.  Thus,  measured  data  and  predictions 
could  be  different  by  a  factor  of  approximately  2  due  to  uncertain 
tv  in  model  meteorological  input  alone. 

C.  DISCUSSION  OF  RESULTS 

This  discussion  is  divided  into  four  sections  --  one  for 
each  of  the  pollutants  measured.  The  included  figures  are  scatter 
plots  of  measured  CO,  NOX,  THC  and  nephelometer  readings  versus 
predicted  concentrat ions .  The  diagonal  lines  drawn  in  these  fig¬ 
ures  enclose  predictions  that  are  within  a  factor  of  two  of  the 
corresponding  measurements.  These  lines  were  found  to  enclose 
greater  than  50$  of  all  the  plotted  points.  Much  of  the  measured 
data  were  invalidated  by  MSI  and  were  therefore  not  available  for 
plotting.  This  is  the  reason  for  the  differences  in  numbers  of 
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plotted  points  from  graph  to  graph  and  was  a  major  limitation 


in  the  validation  effort.  In  fact,  only  approximately  45^  of 
the  pollution  data  taken  during  the  intensive  study  were  con¬ 
sidered  acceptable.  None  of  the  d3ta  taken  during  the  first 
three  days  of  the  intensive  study  met  the  validation  criteria 
in  the  EPA  quality  plan.  However,  those  data  had  to  used  in 
order  to  have  even  a  minimum  of  data  to  compare  with  model  pre¬ 
dictions  . 

Variations  in  predicted  pollutant  concentrations  over  the 
airbase  were  mapped  with  contour  levels  for  the  intensive  study 
and  are  presented  in  Appendix  C.  Contours  for  run  no.  4  (2  Aug, 
1500-1600)  are  included  for  CO  and  PT  concentrations  from  air¬ 
base,  aircraft  and  total  sources.  Contours  for  the  other  nine 
runs  are  included  only  for  CO  and  PT  concentrations  from  air¬ 
craft  sources. 

1 .  Carbon  Monoxide  (CO)  Emissions 

A  comparison  of  the  CO  emitted  each  weekday  with  the 
CO  emitted  during  the  weekend  fperiod  of  reduced  aircraft 
activity)  was  performed  to  better  determine  the  CO  background 
level.  It  was  found  that  on  Saturday  afternoon  the  level  was 
higher  than  that  on  Monday  by  a  factor  of  two,  possibly  due 
to  heaw  traffic  conditions  on  the  surrounding  roadways. 

Also,  on  Sunday,  when  the  winds  were  mostly  calm  or  from  the 
south,  a  high  level  of  CO  was  measured  at  trailer  1.  As  pre¬ 
viously  stated,  weather  conditions  for  the  weekend  during  the 
period  of  intensive  measurement  were  not  representative  of 
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weather  conditions  during  the  weekdays.  Therefore,  no  con¬ 
clusions  could  be  drawn  from  this  comparison  regarding  the 
validity  of  using  trailer  1  measurements  as  indicators  of 
background  CO  levels. 

Figure  7  indicates  that  measured  concentrations  agreed 
with  predicted  total  concentrations  within  a  factor  of  two  at 
trailers  1  and  4.  The  agreement  was  within  a  factor  of  approx¬ 
imately  three  for  trailer  3  data.  (No  measured  CO  data  for 
trailer  2  was  available  during  the  ten  one-hour  time  periods 
used  in  this  study.)  However,  the  good  agreement  may  be  chance 
since  the  environ  input  (land-use  factors,  vehicle  mileage  data, 
etc.)  was  only  estimated.  In  other  words,  what  if  the  high 
levels  of  CO  concentration  at  trailer  1  were  due  to  aircraft, 
but  the  model  did  not  have  properly  input  aircraft  operations 
or  did  not  correctly  determine  dispersion  rates?  AQAM  predicted 
that  the  CO  concentration  due  to  aircraft  at  trailer  1  was 
essentially  zero.  To  check  this,  the  Source  Inventory  program 
was  modified  so  that  all  aircraft  climb  angles  on  takeoff  were 
decreased.  This  maximized  the  near  ground  emissions  from  air¬ 
craft  in  the  area  near  trailer  1.  This  change  had  no  effect  on 
CO  at  trailer  1.  Also,  the  sensitivity  study  discussed  above 
indicated  no  effect  from  increasing  the  hot  refueling  area  and 
hot  refueling  delay  area  source  sizes.  In  other  words,  some 
inaccuracy  in  aircraft  source  specification  near  trailer  1  would 
not  cause  increased  concentrations  at  that  receptor.  Therefore, 
it  appears  to  be  a  valid  assumption  that  trailer  1  was  a  good 
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background  level  indicator  when  a  westerly  wind  prevailed,  and 
the  AQAM  environ  input  for  CO  was  reasonable.  The  model  predicted 
that  CO  concentrations  due  to  environ  sources  were  nearly  constant 
over  the  entire  airbase. 

To  check  the  validity  of  AQAM  predictions  for  CO  emissions  due 
to  aircraft,  trailer  1  measured  concentrations  Cnow  assumed  to  be 
reasonable  background  CO)  were  subtracted  from  the  measured  con¬ 
centrations  at  trailers  3  and  4.  Figure  7  shows  good  agreement  for 
the  very  limited  data  available.  The  higher  predicted  aircraft  CO 
values  at  trailer  3  may  result  either  from  inaccurate  specification 
of  aircraft  idle  CO  emissions  in  the  hot  refueling  area  or  from  a 
too  s lowly- spreading  plume.  It  was  observed  at  MAS,  Miramar  that 
many  times  aircraft  were  operated  at  above  idle  RPM  in  modes  tradi¬ 
tionally  input  into  the  model  as  ''idle”.  To  brieflv  examine  the 
effect  of  the  specified  engine  RPM  on  the  model  predictions,  all 
"idle”  operations  were  changed  to  "normal".  For  most  engines  this 
decreases  the  relative  amounts  of  CO  and  ITHC  while  increasing  MOX 
and  PM.  However,  total  emittants  increase  since  the  fuel  flow  rate 
increases.  In  reality  the  ''normal”  setting  is  far  too  high  to  be 
realistic.  This  high  throttle  setting  is  more  realistic  for  CO 
and  UHC  (since  they  generally  decrease  rapidly  with  RPM  at  low  speeds 
and  then  level  off)  than  for  MOX  and  PT  fwhich  generally  increase 
in  a  linear  manner  with  RPM).  A  more  accurate  method  would  have  been 
to  use  ”h  idle  +  h  normal",  or  etc.  The  increased  power  setting 
improved  the  comparison  with  measured  CO  at  trailer  3  (Fig.  7).  A 
change  of  1  in  stability  category  input  to  AQAM  could  also 
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significantly  change  the  predicted  concentrations  at  trailer  3. 

In  addition  to  predicting  reasonably  accurate  concentra¬ 
tions  at  specific  receptors,  a  model  should  also  correctly 
predict  concentration  profiles  across  the  receptor  grid.  A 
CO  concentration  profile  across  the  airbase  was  constructed 
(Figure  8)  to  illustrate  the  variation  in  predicted  concentra¬ 
tion  along  the  wind  direction.  In  the  two  cases  plotted,  the 
wind  was  from  270°  and  the  stability  category  was  3.  The  two 
profiles  were  plotted  along  the  8  km.  v-coordinate  since  this 
v-coordinate  most  nearly  passed  through  the  trailer  1-4  loca¬ 
tions.  Predicted  and  measured  trailer  data  that  were  available 
were  also  plotted.  "Trailer  profiles"  were  sketched  only  to 
indicate  general  trends  and  do  not  necessarily  represent  actual 
variations.  The  comparison  shows,  as  expected,  that  the  pre¬ 
dicted  trailer  1-4  variation  had  a  much  larger  gradient  than  the 
8  km.  profile  due  to  closer  proximity  to  aircraft  ground  opera¬ 
tions  (taxiways,  hot  refueling  areas,  parking  areas).  The 
measured  profiles  for  both  2  Aug  and  6  Aug  were  similar  to  the 
predicted  profiles,  peaking  between  trailers  2  and  3.  The 
higher  predicted  values  at  trailer  3  were  discussed  above.  Also 
shown  in  Figure  8  is  the  improved  "trailer  profile"  obtained 
when  engine  RPM  was  increased  from  "idle"  to  "normal"  as 
discussed  above. 

2 .  N'OX  Emissions 

Comparison  of  weekend/weekday  data  again  permitted  no 
significant  conclusions  regarding  the  validity  of  using  trailer 
1  as  an  indicator  of  background  NOX . 
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X  COORDINATE  C  KM) 

FIGURE  a 

Predicted  CO  Concentration  Profiles  Along  8  km.  y-coordinate 


Figure  9  presents  measured  versus  predicted  hourly- 
average  NTOX  concentrations  for  trailers  1,  2,  and  4.  (N’o 
measured  data  were  available  for  trailer  5  during  the  ten  one 
hour  time  periods  selected  for  validation  efforts).  As  pre¬ 
viously  stated,  the  comparison  was  based  upon  an  NO-,  conver¬ 
sion  factor  for  ppm  to^gm/m^.  Predicted  concentrations  from 
both  aircraft  sources  alone  and  total  sources  3re  plotted  to 
indicate  their  relative  magnitudes.  Predicted  total  concen¬ 
trations  at  trailers  1  and  4  agreed  with  measured  concentra¬ 
tions  within  a  factor  of  approximately  three.  It  should  be 
noted  that  the  predicted  concentrations  were  all  very  small 
and  varied  much  less  than  the  measured  data.  Also,  the  mea¬ 
sured  data  at  trailer  2  were  much  greater  than  predicted  NOX 
concentrations . 

Because  of  the  general  agreement  between  trailer  1  mea¬ 
sured  and  predicted  concentrations,  it  appears  that  trailer  1 
again  provided  a  good  representation  of  background  concentra¬ 
tions.  Therefore,  trailer  1  measured  concentrations  were 
subtracted  from  those  measured  at  trailers  2  and  4  and  com¬ 
pared  to  predicted  aircraft  NOX  emissions.  Again,  at  trailer 
2  the  measured  (difference)  values  were  much  greater  than  pre 
dieted  aircraft  concentrat ions .  At  trailer  4  the  measured 
(difference)  data  agreed  reasonably  well  with  predicted  air- 
carft  data  (both  were  very  small).  Since  trailer  4  and 
trailer  1  concentrations  were  nearly  the  same  for  both  mea¬ 
sured  and  predicted  data,  and  only  approximately  one-half  of 


the  predicted  trailer  1  values  were  due  to  aircraft,  trailer 
1  was  probably  outside  cost  of  the  aircraft  p luxes  for  the 
existing  wind  conditions. 

Because  trailer  2  was  located  in  a  near-source  region 
where  lateral  concentration  gradients  were  large,  compari¬ 
sons  were  also  made  to  crosswind  receptor  concentrations.  The 
(trailer  2  -  trailer  D!neasurecj  concentrations  were  compared 
to  the  predicted  concentrations  from  aircraft  at  special 
receptor  404  (190m  crosswind/ south  of  trailer  2).  The  pre¬ 
dicted  concentrat ions  were  still  much  less  than  measured 
concentrations,  indicating  that  the  predicted  concentration 
gradients  around  trailer  2  were  not  enough  to  significantly 
improve  the  comparison  between  predictions  and  measurements . 

These  results  indicate  that  the  N'OX  emissions  from  air¬ 
craft  engines  specified  in  AQA.M  are  too  low  for  low  power 
engine  operations  (idle  and  taxi"1.  An  alternative  explana¬ 
tion  is  that  the  aircraft  engine  settings  for  aircraft  located 
around  trailer  2  (hot  refueling  area,  taxiways  ,  and  parking 
areas)  are  well  above  idle,  thus  producing  more  N’OX  than 
assumed  by  AQAM.  Increasing  the  engine  settings  from  idle  to 
normal  as  discussed  above  did  greatly  improve  the  comparison 
between  (trailer  2  -  trailer  1)  predictions  and  measurements 
(from  10  to  136  y^gm/m-5  vs  •  iOO  y^gm/m''  measured  for  2  Aug  from 
1  400-  1  500  hrs o  . 

3 .  Total  Hydrocarbon  (T!^j  Emissions 

The  measured  versus  predicted  total  hourly-averaged 
THC  concentrat ions  for  trailers  1,  3  and  4  are  plotted  in 
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nigure  10.  (Vo  measured  data  were  available  for  trailer  2). 

The  conversion  factor  used  for  ppm  to  /tgm/m0  was  based  on  CH^ 
and  was  therefore  only  an  approximation  for  total  hydrocarbons. 

As  can  ^e  seen  from  the  figure,  predicted  data  were  significant¬ 
ly  lower  and  varied  much  more  than  measured  data.  Measured 
trailer  1  concentrations  were  approximately  1.5  times  greater  than 
trailer  3  concentrations.  This  decrease  is  nearly  the  same  as 
expected  for  downwind  dispersion  from  far  upwind  sources  fi.e., 
due  to  changes  in  C  in  equation  31.  These  results  indicate  that 
almost  all  THC  was  probably  from  environ  sources.  A QAM  predicted 
concentrations  at  trailer  3  were  greater  than  those  at  trailers  1 
and  4  due  to  aircraft  ground  activity.  If  most  of  the  measured 
concentrations  of  THC  are  in  fact  due  to  environ  sources  and 
measured  trailer  1  values  are  accurate,  then  either  AQAM  values 
for  THC  emittants  due  to  environ  sources  are  low  (i.e.,  land-use 
factors  are  low")  or  the  conversion  factor  toAgm/m0  is  in  signifi¬ 
cant  error.  The  former  would  also  imply  that  the  values  used  in 
AQAM  for  THC  emittants  from  aircraft  sources  are  too  high  (at 
trailer  3  downwind  of  the  hot  refueling  area''.  This  particular 
observation  could  have  been  better  clarified  had  measured  data 
been  available  from  trailer  2.  Increasing  engine  RPM  as  discussed 
above  in  this  case  reduces  THC  from  aircraft.  When  combined  with 
increased  environ  sources  of  THC,  this  change  would  improve  the 
agreement  with  measurements. 

4 .  Particulate  (PT1  emissions 

Figure  11  is  a  plot  of  the  measured  (converted  bscat' 
versus  predicted  total  hourly  averaged  PT  concentrations.  The 


measured  data  were  within  -  40",  of  the  mean  value.  The  mea¬ 
sured  values  at  trailers  1  and  4  were  essentially  the  same. 

The  comparison  is  fairly  good  (within  a  factor  of  three  for 
“O',  of  the  data)  at  trailers  1  and  4  using  the  aforementioned 
conversion  factor  for  bscat  to  ^gm/m0.  The  model,  however, 
appears  to  overpredict  PT  concentrations  at  trailer  2.  A.QAM 
predicts  that  most  of  the  PT  concentration  is  from  aircraft 
sources.  Therefore,  if  trailer  1  data  are  good  indicators  of 
background  PT  concentration,  then  AQAM  has  low  environ  source 
PT  input  fland-use  factors,  vehicle  mileage,  etc.)  and/or  high 
aircraft  source  PT  input.  Increasing  engine  RPM  as  discussed 
above  increases  PT  and  makes  the  comparison  generally  less 
favorable. 

C.  CONCLUSIONS  AND  RECOMMENDATIONS 

Approximately  SO"  of  the  predicted  levels  of  concentra¬ 
tion  were  found  to  agree  with  measured  levels  within  a  factor 
of  two.  The  results  also  indicated  that:  (11  predicted  CO 
concentrations  agreed  quite  well  with  measured  data;  (2)  model 
predictions  were  too  low  for  NOX  emissions  from  aircraft  oper¬ 
ating  in  the  idle/taxi  mode;  and  (3)  predicted  THC  and  PT  con¬ 
centrations  were  too  high  for  aircraft  operating  in  the  idle/ 
taxi  mode  and/or  were  too  low  for  environ  sources.  The  latter 
appears  more  reasonable.  Agreement  between  model  predictions 
and  measured  values  was  significantly  improved  by  increasing 
engine  RPM  settings  above  idle  in  all  modes  normally  specified 


6  3 


at  idle. 


For  a  reasonably  complete  model  validation  to  be  accom¬ 
plished  much  more  measured  data  must  be  obtained  during  a 
specific  time  period  of  observed  meteorological  and  operation¬ 
al  activity.  The  conclusions  from  this  intensive  study  were 
based  on  very  limited  data  and  can  only  be  considered  prelim¬ 
inary  results.  Accurate  data  for  background  levels/local  air 
quality  are  important  for  determination  of  aircraft  source  con¬ 
tributions  to  total  emittants.  It  would  be  most  beneficial  to 
obtain  pollution  measurements  on  weekends  at  a  time  when  air¬ 
craft  activity  is  low  and  meteorological  conditions  are  very 
similar  to  weekday  conditions.  If  at  all  possible,  any  additional 
intensive  efforts  should  be  conducted  during  a  period  with  less 
variations  in  meteorology.  Detailed  data  collection  should  begin 
several  days  before  the  detailed  operational  data  are  collected 
in  order  to  ensure  a  more  complete  data  set  than  was  obtained  in 
this  initial  effort. 

Model  predictions  were  very  sensitive  to  the  specified 
stability  category.  Model  validation  efforts  could  be  improved 
if  stability  classes  could  be  measured  and  specified  to  half¬ 
integer  values. 

Jet  dispersion  rate  differences  from  those  for  elevated, 


low  velocity  sources  and  variations  with  orientation  to  the 
wind  direction  require  further  study.  Plume  rise  of  iet  exhausts 
should  also  receive  additional  study. 
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AIRCRAFT  FT  CONCENTRATION  PROFILE  (RUN  NO.  1 ) 
(Scale  =  20  per  contour) 


AIR3A3E  CO  CONCENTRATION  PROFILE  (RUN  NO.  1  ) 
(Scale  =  1  //gm/m^  per  contour) 
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OTAL  CO  CONCENTRATION  PROFILE  (RUN  NO.  1  ) 
(Scale  =  50  per  contou. ) 
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TOTAL  PT  CONCENTRATION  PROFILE  (RUN  NO.  1  ) 
(3cale  =  20  >»gra/m^  per  contour) 
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AIRCRAFT  PT  CONCENTRATION  PROFILE  (RUN  NO.  2) 
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AIRCRAFT  CO  CONCENTRATION  PROFILE  (RUN  NO.  3) 
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AIRCRAFT  CO  CONCENTRATION  PROFILE  (RUN  NO.  4) 
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AIRCRAFT  CO  CONCENTRATION  PROFILE  (RUN  NO.  6) 
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AIRCRAFT  PT  CONCENTRATION  PROFILE  (RUN  NO.  7) 
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AIRCRAFT  PT  CONCERT RAT ION  PROFILE  (RUN  NO.  S) 
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AIRCRAFT  CO  CONCENTRATION  PROFILE  (RUN  NO.  9) 


AIRCRAFT  PT  CONCENTRATION  PROFILE  (RUN  NO.  9) 


APPENDIX  C 

CO  AND  ?T  CONCENTRATION  PROFILES  FROM  AIRCRAFT  SOURCES 

(INTENSIVE  STUDY) 
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AIRCRAFT  FT  CONCENTRATION  PROFILE  (1  AUG  1  300-1400) 
INCREMENTED  FROM  30.0 
(Scale  =  30  ^ugm/m^  per  contour) 


AIRCRAFT  CO  CONCENTRATION  PROFILE  (1  AUG  1L00-1500) 
INCREMENTED  FROM  50.0 
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AIRCRAFT  PT  CONCENTRATION  PROFILE  (1  AUG  1 UOO-1 500) 
INCREMENTED  FROM  30.0 

ini 


AIRCRAFT  PT  CONCENTRATION  PROFILE  (2  AUG  II4.OO-15OO) 
INCREMENTED  FROM  30.0 
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AIRCRAFT  CO  CONCENTRATION  PROFILE  (2  AUG  1 500-1 600) 
INCREMENTED  FROM  50.0 


AIRBASE  GO  C  ON  C  SOT  SAT  I C II  PROFILE  (2  AUG  1  500-1  600) 
INCREMEJITED  FROM  1.0 
(Scale  =  U  per  contour) 


AIRBASE  PT  CONCENTRATION  PROFILE  (2  AUG  1 500-1 600) 
INCREMENTED  FROM  10.0 
(Scale  *  5  per  contour) 


J 


106 


TOTAL  PT  CONCENTRATION  PROFILE  (2  AUG  1 500-1 600) 
INCREMENTED  FROM  50.0 
(Scale  =  50  jj gm/m^  per  contour) 


AIRCRAFT  CC  CONCENTRATION  PROFILE  (2  AUG  1515-1615 
INCREMENTED  FROM  50.0 


AIRCRAFT  PT  CONCENTRATION  PROFILE  (2  AUG  1515-1615) 
INCREMENTED  FROM  30.0 
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AIRCRAFT  CO  CONCENTRATION  PROFILE  (3  AUG  11 00-1 200) 
INCREMENTED  FROM  50.0 
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AIRCRAFT  ?T  CONCENTRATION  PROFILE  (6  AUG  II4.OO-1 500) 
INCREMENTED  FROM  30. 0 
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AIRCRAFT  PT  CONCENTRATION  PROFILE  (6  AUG  1500-1600) 
INCREMENTED  FROM  30.0 


AIRCRAFT  CO  CONCENTRATION  PROFILE  (6  AUG  1515-1 61 5) 
INCREMENTED  FROM  50.0 
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AIRCRAFT  PT  CONCENTRATION  PROFILE  (6  AUG  1515-1615) 
INCREMENTED  FROM  30.0 
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AIRCRAFT  CO  CONCENTRATION  PROFILE  (7  AUG  1 500-1 600 ) 
INCREMENTED  FROM  50.0 
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AIRCRAFT  PT  CONCENTRATION  PROFILE  (7  AUG  1 500-1 600) 
INCREMENTED  FROM  30.0 
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